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We study front propagation in stirred media using a simpli- 
fied modelization of the turbulent flow. Computer simulations 
reveal the existence of the two limiting propagation modes 
observed in recent experiments with liquid phase isothermal 
reactions. These two modes respectively correspond to a wrin- 
kled although sharp propagating interface and to a broadened 
one. Specific laws relative to the enhancement of the front ve- 
locity in each regime are confirmed by our simulations. 

PACS numbers: 05.40+j, 47.27.Sd 



I. INTRODUCTION 

Front propagation problems have been receiving an 
increasing attention due to its relevance for non- 
equilibrium studies in a wide variety of fields [Q. Most 
of these studies are concerned with the velocity selec- 
tion problem in deterministic systems or homogeneous 
medium (reaction-diffusion equations) . Very recently 
front propagation in modulated |^ or stochastic media 
has been discussed Q]. 

Self sustained chemical fronts propagating in quies- 
cent media advance steadily at a velocity vq which re- 
sults purely from the interplay between the chemical 
time scale Tchem and the species diffusivity coefficient 
D: Vq ~ {^'''chemy^^ ■ Howcvcr, in stirred media, addi- 
tional length and time scales come to play a fundamental 
role into such a dynamical process. The interface, either 
wrinkled or broadened by the turbulent flow, propagates 
at a velocity vt which is enhanced relative to vq to an 
extent which depends not only on the invested energy 
but also on the coupling between the intrinsic length and 
time scales and those externally introduced by the stir- 
ring mechanism [^. 

Needless to say, the addressed problem is of singular 
importance either as a way to inquiry about fundamental 
aspects of turbulent flows, but mainly for its practical rel- 
evance to combustion processes However, abundant 
studies in either way motivated have not been able yet 



to resolve the existing considerable discrepancies between 
theoretical predictions and experimental results. These 
discrepancies concern fundamental issues such a turbu- 
lent propagation rates, velocity quenching effects, role of 
turbulent spectra, etc. §. 

From the experimental point of view much has been 
gained recently by enlarging the scope of the problem 
when going from combustion processes to liquid phase re- 
actions . In this way experimentalists more confidently 
approach the set of simplifying assumptions invoked by 
most of the theoretical models. These later ones, either 
based on the so-called G-equation H , passage-rate anal- 
ysis or scale-covariance laws [^0[ tend to focus on a 
single mode of turbulent propagation i.e. that identi- 
fied later on as based on the so-called mechanism. On 
the other hand, when resorting to computing simulations 
pd| , numerical accuracy commonly restricts the exam- 
ined scenario to small stirring intensities. On the whole, 
it is then likely to expect a future demand for more versa- 
tile computer approaches to compare with more general 
and better controlled experimental observations. 

Advancing towards that direction, we report here on 
our computer results for interfaces propagating in stirred 
media. Stirring is mimicked by using an stochastic differ- 
ential equation (SDE) based algorithm to simulate ran- 
dom, stationary, isotropic and statistically homogeneous 
fiows with well prescribed properties ||l2| . In spite of this 
"ad hoc" procedure, the pair of experimentally identified 
limiting modes of front propagation in liquids phases , 
i.e. the aforementioned Huygens propagation (HP) and 
the distributed reaction zone (DRZ) regimes, are clearly 
exhibited in our simulations. In addition we do not suf- 
fer from strict limitations concerning stirring intensity, 
what enables us to look for broad range scalings of the 
turbulent propagation rates. Finally, the main advan- 
tage of our model is the easy control we gain over the 
fiow parameters more relevant for front propagation in 
stirred media i.e., the intensity Uq, spectra E{k), spatial 
correlation length Iq, and time correlation t^. In partic- 
ular on what follows we adopt the Kraichnan spectrum 
]l3[ , describing a widely distributed band of excitations 
around a peak centered at some well-defined wavenumber 
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ko At this point is worth emphasizing that the proposed 
methodology is totally adapted to a statistical descrip- 
tion of turbulence (l4j and belongs to the class of the 
so-called synthetic turbulence generating models jisl] . 

The organization of this paper is as follows. The 
reaction-diffusion-advection model and the existing ana- 
lytical prediction on front propagation in turbulent media 
are briefly reviewed in Sec. Our numerical results are 
presented and commented according the above theoreti- 
cal predictions in Sec. 



II. TURBULENT FRONT MODEL AND 
ANALYTICAL PREDICTIONS 

A detailed presentation of the algorithm can be found 
in , but it is worth summarizing here its main ingre- 
dients: 

i) A stream function ri{r, t) is introduced trough a SDE 
written in terms of a zero mean Gaussian white noise 



dt 



Q[A2v2]V.C(r,i), 



(1) 



(C(ri,ii)C^'(r2,i2)) = 2e^v5{t^-h)5{v^-V2)5''. (2) 

ii)From the stochastic field 77, the 2-d incompressible 
advecting flow is obtained as usual 



v(r,i) 



drj{v,t) dri{r,t) 



dy 



dx 



(3) 



The two dimensional space is chosen here for several rea- 
sons. First of all, it is chosen for the sake of simplicity 
of numerical simulations. Secondly because the essential 
trends of the phenomenology we want to study are inde- 
pendent of the spatial dimension. Recent experiments in 
chemical reactions in quasi 2-d geometries [Q and theo- 
retical analysis [ p^ support this assumption. 

iii) The whole procedure is Fourier transformed and 
accordingly discretized. 

iv) Within this scheme, the homogeneous, isotropic 
and stationary velocity correlation R{r, s) and the en- 
ergy spectrum E(k) can be properly defined and easily 
evaluated jl^. From these quantities, the stirring inten- 
sity Mg, and the integral time to and length Iq scales of 
the flow can be obtained using the following standard 
prescription: 



(4) 
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The adopted Kraichnan spectrum has the following 
form 

cx fc^exp [-A^fc2] . (5) 
This energy distribution is reproduced by prescribing 



Q(A^V^) = cxp 



(6) 



On the other hand, the basic statistical parameters of 
the flow introduced above expressed in terms of the input 
parameters, eo,A, z^, read 



SttA^ v 2 



(7) 



The next step is the formulation of the reaction- 
diffusion-advection scheme. We model this situation by 
means of a dynamical equation for the passive scalar field 
'(/'(r, t) in a two dimensional space: 



V.(vV'). 



(8) 



where f{tp) is a nonlinear reaction term with at least two 
steady states. In our simulation we have taken /(ip) = 
ijj^ — Our results are independent of the particular 
form of f{tp) chosen in the simulations. 

In the absence of stirring (v = 0), a stable planar front 
propagates the stable state tp = 1 ("products") into the 
invaded unstable one ip = ( "reactants" ) . The dimen- 
sionless front velocity and front thickness respectively 
read: 




(9) 



Superposing stirring, two limiting regimes of front 
propagation are found in our simulations. First when 
the typical length scale of the flow Iq is larger than the 
intrinsic one associated to the reaction-diffusion dynam- 
ics So, we observe a distorted front which propagates as 
a still rather sharp interface. Such a propagation mech- 
anism is known in the combustion literature as the "thin 
flame," "flamelet" or "reaction sheet" regime 0. Con- 
trarily when Iq is smaller than So, we observe what is 
referred in the literature as a "DRZ" regime i.e., a 
broadened front disrupted by the stirring flow. In both 
cases, turbulent propagation rates are larger than in qui- 
escent media (parameter values of the presented simula- 
tions are quoted in the figure captions). 

Each one of the previously identified conditions corre- 
spond to a specific propagation mechanism. The com- 
mon rationale behind the "thin flame" mode is based on 
a HP-like argument: the front has the same local struc- 
ture as in the planar case with normal velocity given by 
Wo, but its length increases due to wrinkling. This results 
on faster propagation velocities, in such a way that the 
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relative increment in the velocity is equal to the relative 
increment of the length. On the other hand stirring is as- 
sumed to affect the velocity in the DRZ regime by solely 
increasing diffusive transport inside the broadened front. 
Note in passing that in neither case any effect of the stir- 
ring on the intrinsic chemical time scale is considered. 

Before to translate these arguments into quantitative 
terms let us define precisely the quantities involved. Lq 
will be the length of the planar front which coincides 
with the lateral length of the system. Lt is the length of 
the curved front evaluated numerically as the length of 
the curve level of the field at the value = 1/2. Dt is 
the effective diffusion when turbulent flow is present, ft 
is numerically evaluated by simulating the diffusion of a 
passive scalar under the influence of the same flow free of 
any chemical front. The velocity of the front under the 
turbulent flow vt is evaluated numerically as the incre- 
ment of the products per unit of time and normalized to 
Lo- 

According to our previous arguments we can establish 
the following two analytical results: 
i) In the HP regime, 



Urp 

Vo ~ 



S 



(10) 



ii) Contrarily, in the DRZ mode, we simply adapt the 
first fundamental relation of Eq. (O) to obtain 



Dt 
~D 



(11) 



III. NUMERICAL RESULTS AND COMMENTS 

The next and most involved step consists, however, in 
using Eqs. ( p^ ) and (|l^) to compare vt, or its dinien- 
sionless form S, as a function of the stirring intensity 



or its dimensionless value 



Q with the results 



u-o, 

obtained from the numerical simulations of the Eq. (8). 
Let's consider separately on what follows the HP and 
DRZ regimes. 



A. HP mode 

On what refers to the HP mode, our first task was 
to check relation ([lO|). The collected data for the differ- 
ent values of Wq are summarized in Fig. 1. For the sake 
of comparison, we include in this figure results obtained 
from two additional and somewhat related stirring condi- 
tions. The first one, hereafter referred as frozen stirring, 
corresponds to a fixed configuration of the random flow. 
The second one, referred on what follows as periodic stir- 
ring, represents nothing but the limit of a deterministic 



and single scale flow, constructed from the single mode 
stream function 



n-KX niry 
y) = m cos(-^) cos(-^). 



(12) 



representing a periodic array of n x n eddies, where L 
is the system size. According to Fig. 1, the geometric 
argument leading to (|l^) seems well-supported by our 
simulations with the simple exception of those situations 
involving very intense periodic flows, were the presence of 
overhangs (leading to the formation of isolated islands of 
reactants) is unavoidable. Note in this respect that the 
time evolution of those islands contribute positively to 
the computed velocity, measured as the time variation of 
the rate of occupation of the -0 = 1 state, but negatively 
to the front length. 



B. DRZ mode 

Results for the DRZ propagating mode are summa- 
rized in Fig. 2. Here again results for 5* obtained by 
directly simulating the front propagation dynamics un- 
der the three different stirring mechanisms so far con- 
sidered are plotted together. Computer results are com- 
pared with the corresponding theoretically predicted val- 
ues based on Eq. (|ll|) above. The agreement based on 
Eq. (^l|) is remarkable for the whole range of Q values 
here considered and irrespective of the stirring flow. Ex- 
plicit theoretical dependences of Dt on Uq valid respec- 
tively in two asymptotic cases are: 

i) For a random flow in the weak stirring limit and 
small integral time to 



Dt- D 



(13) 



A result obtained from Eq. (8) considering v as a Gaus- 
sian white noise 

11] and, 

ii) For a periodic flows in the limit of small Peclet num- 
ber 0, 

Dt - D ^ ul/D. (14) 

These two analytical results used in relation to Eq. (|Tl| ) 
above lead to the theoretical estimates also shown in 
Fig. 2 (continuous lines). 

A last remark concerning these two flgures can clar- 
ify these results. The argument leading to Eq. (10) and 
Fig. 1 is pure kinematical or geometrical and only de- 
pends on the fact that the front is narrow enough. Thus 
different dynamical situations flt the prediction once this 
condition is fulfilled. Contrarily the argument for the 
DRZ (Fig. 2) case is dynamical so we observe different 
behaviors for the three different flows. 

In summary, a recently proposed algorithm to model 
random stirring has been used in relation with front prop- 
agation. Although artiflcial in nature such a model re- 
produces realistic scenarios of front propagation modes. 
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These conchisions together with the versatile computer 
implementation of the proposed algorithm give us new 
hopes to address some of the open issues in the field of 
turbulent front propagation. Research in this direction 
basically aimed at finding the power law dependences of 
the propagating velocities on the stirring intensities and 
their fits to some of the most well known relations pro- 
posed in the literature, are being presently conducted and 
will be reported elsewhere. 
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FIG. 1. S versus Lt/Lq (Eq. 10) for turbulent flow (cir- 
cles, lo = 4.0, to = 1.0), random frozen (squares, lo = 4.0) 
and periodic eddies (triangles, 8x8). In all the simulations 
we have employed a square lattice of 128 x 128 points and 
unit spacing Ar = 0.5 and D = 0.3. 
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FIG. 2. S versus for the DRZ regime. Full symbols are 
simulation results of the front propagation while open symbols 
stand for the theoretical Eq. (11) using the effective turbulent 

diffusion (sec text). Continuous lines correspond to the an- 
alytical predictions based on the perturbative results in (13) 
and (14). Circles, squares, and triangles correspond to turbu- 
lent flow {lo — 2.0, to = 0.1), random frozen {Iq — 2.0) flows 
and periodic eddies respectively (16 x 16 eddies) {D = 2.0). 
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